Evidence for granulation in early A-type stars. 
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ABSTRACT 

Stars with spectral types earlier than about FO on (or close) to the main sequence have long 
been believed to lack observable surface convection, although evolutionary models of A-type stars 
do predict very thin surface convective zones. We present evidence for granulation in two 5 Scuti 
stars of spectral type A2: HD 174936 and HD 50844. Recent analyses of space-based CoRoT 
(Convection, Rotation, and planetary Transits) data revealed up to some 1000 frequencies in the 
photometry of these stars. The frequencies were interpreted as individual pulsation modes. If 
true, there must be large numbers of nonradial modes of very high degree I which should suf- 
fer cancellation effects in disk- integrated photometry (even of high space-based precision) . The 
p-mode interpretation of all the frequencies in HD 174936 and HD 50844 depends on the assump- 
tion of white (frequency independent) noise. Our independent analyses of the data provide an 
alternative explanation: most of the peaks in the Fourier spectra are the signature of non-white 
granulation background noise, and less than about 100 of the frequencies are actual stellar p- 
modes in each star. We find granulation time scales which are consistent with scaling relations 
that describe cooler stars with known surface convection. If the granulation interpretation is 
correct, the hundreds of low-amplitude Fourier peaks reported in recent studies are falsely in- 
terpreted as independent pulsation modes and a significantly lower number of frequencies are 
associated with pulsation, consistent with only modes of low degree. 

Subject headings: stars: oscillations stars: individual: HD 50844, HD 174936 stars: variables: delta 
Scuti 



Do more frequencies actually mean 
more modes? 



The (5 Scuti stars are intermediate-mass pulsat- 
ing variables with spectral types A — F in the 
lower part of the classical instability strip. They 
are opacity-driven pulsators, which can exhibit 
large numbers of excited acoustic modes, mak- 
ing them key candidates for asteroseismology. Be- 
cause of this modeling potential, huge efforts have 
been made to observe 5 Scuti stars from ground, 
offering increasingly rich eigenfrequency spectra, 
thanks to steady improvements of observational 
techniques. The pinnacle o f this attemp t s was 
the star FG Vir, for which iBreger et all (|2005l ) 
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detected 75 frequencies with amplitudes down to 
^200 ppm. The trend of finding oscillations of 
lower and lower amplitude led to the belief that, 
only the tip of the eigenspectrum iceberg was visi- 
ble and as noise levels were pushed lower by better 
observations, more modes would be found in most 
pulsating stars. 

This has been impressively confirmed by re- 
cent space-based observations. iMatthewsl (|2007j) 
reported on MOST (Micrqvariability fc Oscilla- 



tion of STars; Walker ct al. 2003; Matthews ct al 



2004) observations of HD 209775, revealing a sim- 



ilar number of frequencies as had been found in 
FG Vir. But truly groundbreaki ng results came 



from CoRoT ([Baglin et al.l 120061 ) , which uncov 



ered an unexpectedly rich eigenspectrum in the 
S Scuti star HD 50 844, with potentially more than 
1000 



i star tiL) 0U 844, wit n potentially more tnan 
frequencies ( Poretti et al.l 20091) . Similar 



was f ound for HD 174936 (jGarcia Hernandez et al 
2009), and results on other 5 Scuti stars in the 



CoRoT sample are currently processed. The ba- 
sic result is that several hundreds of frequencies 
have been detected in several of the S Scuti stars 
monitored with CoRoT. 

This early result inspired a new perspective 
on 5 Scuti pulsation, potentially addressing long- 
standing problems of the mode amplitude distri- 
bution and unknown mode selection process(es) 
in 6 Scuti stars. But at the same time, it also 
opened new fundamental problems: how can such 
large numbers of modes be actually detected in 
measurements of integrated light from stars? As- 
teroseismic models predict some tens of acoustic 
modes (maybe one hundred if rotationally split 
frequencies are included in the tally) with de- 
grees I < 4, which are potentially excited in the 
frequency range of S Scuti-type oscillations. A 
substantially larger number of modes than this 
would have to include modes of degree I > 4 
whose net amplitudes should suffer strong can- 
cellation effects wh en integrated across the stel- 



lar d isk (see, e.g., iDaszvhska-Daszkiewicz et al 
20061 ). We find it implausible that hundreds of 



modes with high degree I would survive amplitude 
cancellation, even with space-based precision, un- 
less some effective spatial filter was present on the 
stellar disk. High-contrast large-scale structures 
on the stellar surface could act as such a spacial 
filter. Such features are known for active cool stars 
(and for strongly magnetic B — F stars) but have 
not yet been observed among S Scuti stars. We 
prefer a simpler interpretation for the multitude 
of frequencies in the newly observed power spec- 
tra of 5 Scuti stars: non- white stellar background 
noise, due to surface granulation. 

Stars in the 8 Scuti part of the HRD are believed 
to have no measurable surface convection. Be- 
cause of this long-standing view, space-based pho- 
tometry of 5 Scuti stars has been analysed assum- 
ing the noise is dominated by white (frequency- 
independent) photon noise. The CoRoT light 
curves of HD 174936 and HD 50844 were processed 
using standard iterative sine-wa ve fitting algo - 



mg 

rithms. One of them, SigSpec (jReegenl |2007( ). 
assumes the variance in the time series to be 
dominated by white noise and rates the signifi- 
cance of the signal components according to this 
assumption. This approach clearly fails in the 
case of frequency-dependent noise. The noise 
peaks in the frequency range of an enhanced noise 



can have much larger amplitudes than the real 
(white) noise and are potentially misinterpreted 
as significant signal. Other tools, like Period04 
(|Lenz fc Bregerll2005l ). are less restrictive and can 
handle frequency-dependent noise. On the other 
hand, the program has technical difficulties to si- 
multaneously fit several hundred signal compo- 
nents and is therefore used to check/confirm the 
main frequencies (up to say 100 in number) iden- 
tified by other programs. 

Non-white intrinsic stellar noise is a feature of 
cool stars with convective envelopes. The tur- 
bulent motions in the stars' outer convective en- 
velopes generate quasi-stochastic power acting on 
various time scales, with amplitudes strongly de- 
creasing with increasing frequency. Although the 
signal is stochastic, it exhibits distinctive charac- 
teristics. The Fourier transform of the signal fol- 
lows simple power laws where different physical 
processes on (or near) the stellar surface produce 
the same type of signal but on different time and 
amplitude scales. For the Sun, the different signal 
components are usually assigned to stellar activity, 
activity of the photospheric/chromospheric mag- 
netic network, and granulation, with time scales 
ranging from months (for act ive regions) to min - 
utes (for granulation) (see e.g. . [Michel et al.l 2009). 
This has also been seen i n other Sun-like stars 
(see e.g. JMichel et al . 2008), and also for the more 
evolved red giants (see e.g., Kallinger et al.l l2009). 

The presence of a granulation background sig- 
nal depends of course on the presence of a sur- 
face convection zone — a property which is usu- 
ally attributed to cool stars beyond the red bor- 
der of the instability strip. But according to stel- 
lar evolutionary models, stars located within the 
instability strip have a thin, but nevertheless non- 
negligible, convective surface layer (see Sect. 3). 
The proper modeling of convection is still a weak 
facet of most stellar evolution calculations. But 
the existing models give reason to consider that 
a subtle granulation background signal should be 
observable in A stars. We believe that this was not 
observed (at least to our knowledge) and recog- 
nised before because the amplitudes are very low, 
requiring precise, long, and uninterrupted obser- 
vations to detect the signal. Granulation ampli- 
tudes are measured in power density because the 
actual amplitude in the observations depend on 
the sampling and length of the time series. For 
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the Sun, granulation has an a mplitude on the or 



der of less than lppm 2 / /xHz (| Michel et al. 
and typically 10 3 to 10 4 p pm 2 / /xHz for stars o n 



iooi) 



the lower red-giant branch ( Kallinger et al. 2009). 
We expect the amplitudes for (5 Scuti stars to be 
somewhere in between. 



Power spectrum modeling: 
background and pulsations 



Stellar 



For the Sun, it is common practice to model 
the background signal with a sum of power laws, 
P{ v ) = Si a */(1 + { v /bi) Ci ), where v is the fre- 
quency, aj, hi, and Cj are the amplitude, charac- 
teristic frequency (inverse time scale), and slope 
of the power laws, respectively. The number of 
components usually ranges between two and five, 
depending on the frequency coverage of the obser- 
vat ions. Power-l aw models were first introduced 
by Harvcv (1985) with t he slope fixed to 2. But 



Aigrain et ail (|2004h and lMichel et al.l (|2009t ) have 
shown that, at least for the Sun, the true slope is 
closer to 4. The presence of an additional signal, 
like solar-type p modes, will significantly distort 
the fit, and one has to include a ter m in a model to 
accou nt for that additional signal. iKallinger et al 



(2009) lave approximated the shape of the pul- 
sation power excess in red giants by a Gaussian, 
which is also r easonable for main-sequen ce solar- 
type pulsators ( Gruberbauer et aTl 12009). This is 
more complicated for 5 Scuti stars as the excess 
envelope can be quite different from a Gaussian. 
We therefore pre- whiten the largest amplitude fre- 
quencies and assume the residual power excess to 
be at least Gaussian-like. For a detailed frequency 
analysis a more realistic model might be needed. 
But this is not our current goal. 

Following IKallinger et al. I (<2009h . we model the 
residual power spectra with a superposition of 
white noise, a sum of power laws, and a Gaussian. 



(1) 

where P n represents the white noise contribu- 
and <j g are the height, cen- 



tion and P g , 
tral frequency, and width of the power excess 
hump, respectively. We use a Bayesian Markov- 
Chain Monte Carlo algorithm to fit the global 



model to the residual power spectra. Details 
about the fitting proce dure may be found in 



Gruberb auer et al. (2009) 
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2.1. HD 174936 

The field 5 Scuti star HD 174936 (A2; V = 
8.6) was observed within the CoRoT asteroseis- 
mology program during the first 27-day short 
ru n campaign. A first analysis was presented 
by iGarcia Hernandez et al" ( 20091) . who detected 
about 420 individual frequencies with SNR > 4.0, 
assuming white noise. This corresponds to a de- 
tection limit of about 14 ppm. 

We use SigSpec to pre-whiten the ten most 
significant frequencies from the CoRoT time se- 
ries, which are identic al to t hose reported by 
Garc ia Herna ndez et al Although the 

star was monitored continuously, a significant 
number of data points are flagged d uring the re- 
duction process ( Samadi et al.l 20071) . mainly out- 
liers due to high-energy particle impacts during 
the satellite's pa sses through the South A tlanti c 
Anomaly. Like Garcia Hernandez et al.l (2009), 
we simply remove these questionable data points, 
degrading the duty cycle of the time series to 
^90%. The resulting gaps are fairly regular, 
producing low-amplitude aliasing peaks in the 
power spectrum. This is not a problem for a 
pre-whitening sequence but distorts a direct fit to 
the spectrum. We therefore fill the gaps by linear 
interpolation, which is the sta ndard procedure for 



solar -type pulsators (see e.g. lAppourchaux et al 
2008). The power spectra of the original data 



and of the prc-whitencd and gap-filled data are 
plotted in Fig.[TJ with the spectral window of the 
original data shown in the inset. The gap-filling 
suppresses any alias signal and obviously changes 
the power spectrum dramatically for frequencies 
above ~600 //Hz. Note that such gap-filling might 
influence the determinations of individual frequen- 
cies and amplitudes, which is not our goal in this 
study. 

We then fit a global model (Eq.[2|) to the resid- 
ual power spectrum. In the case of HD 174936, 
the residual S Scuti-type oscillation power appears 
to be concentrated in a single power excess cen- 
tered on about 350 /iHz. The resulting global fit 
reproduces the overall shape of the residual power 
spectrum and demonstrates the strong frequency 
dependence of the background signal. At low fre- 
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Frequency ((iHz) 

Fig. 1. — top panels: Original (grey curve) and residual gap-filled (black curve) spectra of HD 50844 (right) 
and HD 174936 (left) after pre-whitening the 20 and 10 most significant frequencies, respectively, and a global 
model fit (blue curve - in the online version only). Red lines represent the models without the pulsational 
components consisting of two power laws (dashed lines) and white noise, bottom panels: Residual power 
spectra normalised to the power-law components. Dashed lines indicate a signal-to-background-noise ratio 
of SNbR = 4.0. 



quencies, we find the background signal to be as 
high as about 36ppm. It decreases to ^19ppm 
at 100 /xHz before it fades into the white noise of 
^2 ppm. 

Our next step is to use the power-law compo- 
nents of the fit to correct the power spectrum for 
the frequency-dependent background noise. The 
resulting signal-to-background-noise ratio (SNbR) 
is given in Fig.[T] We find 55 frequencies exceed- 
ing a significance limit of SNbR > 4.0 resulting 
in a total of 65 significant frequencies, which we 
attribute to d Scuti-type pulsations. This is small 
co mpared to the 420 pulsa t ion m odes identified 
by iGarcfa Hernandez et al. ( 2009f ). but in much 
better agreement with the maximum number of 
frequencies reasonably expected for low-degree p- 
modes in a S Scuti star. 

HD 174936 is a metal-deficient star with an 



iron a bundance [Fe/H] ~ —0.32 (jCharpinet et al 



2006). We have determined an effective tempera- 
ture 7^g_^_8080±160 K from Str omgren photome- 
try (Irlauck et al.lll998 ) using the Moon fc Dworetskv 
(|1985l ) calibration, which gives also an absolute vi- 



sual magnitude My = 1.91±0.32. With M bol Q 
= 4.75 and a bolometric c orrection BCy — 0.003 
( Lejeune fc Schaererl 12001 ). we obtain L/LQ = 
14.4 ± 4.4 and R/R® = 1.91±0.35. 

2.2. HD 50844 

HD 50844 (A2; V = 9.09) was continuously ob- 
served by CoRoT during the 58-day initial run. 
A detailed frequency analysis and fi rst asteroseis- 
mic in terpretation was performed bv lPoretti et al 



(|2009h . In their analysis, they must pre- whiten 
the time series with more than 1000 individual 
frequency terms to decrease the residual rms in 
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the time series to what they believe is the detec- 
tion limit of the data. From their Fig. 3, it is ob- 
vious that they would need yet thousands more 
frequencies to bring the mean amplitude at fre- 
quencies below 50 to the mean amplitude at high 
frequencies. This alone is a strong argument for 
the presence of strongly non- white noise. 

Following the approach described above, we 
pre-whitcn the 20 most significant frequencies 
from the CoRoT data of HD 50844 and fill the 
gaps in the residual time series by linear interpo- 
lation. Unlike HD 174936, the oscillation power of 
HD 50844 seem to be distributed over a broader 
frequency range, with two clusters around about 
50 and 180 fiHz. We therefore include a second 
Gaussian in our global model of the power spec- 
trum (Eq.[2]). The power spectrum of the origi- 
nal data is shown in Fig.Q] along with the power 
spectrum of the pre-whitened and gap-filled data, 
as well as the most probable global model with 
and without the Gaussian terms. In the case of 
HD 50844, the contrast between the low-frequency 
background signal and the high-frequency white 
noise is even higher than for HD 174936. We find a 
low-frequency background level of about 160 ppm, 
which decreases to ~35ppm at 100 /iHz, before 
it fades into the high-frequency white noise of 
^2 ppm. (Coincidentally, the white noise com- 
ponent for HD 50844 is almost the same as for 
HD 174936. This is because the average photon 
flux of the HD 50844 observations is smaller than 
for HD 174936, but the time series consists of more 
than twice as many data points.) 

We find 105 frequencies exceeding a SNbR > 
4.0, resulting in a total of 125 significant fre- 
quencies which we attribute to 6 Scuti-type pul- 
sat ions, compar e d to more than 1000 suggested 
bv lPoretti etaH (|2009h . 

HD 50844 is also a metal-poor star, slightly 
more evolved than HD 174936, with T fiff = 7500±20 0K, 
M v = 1.31, and [Fe/H] ~ -0.4 dPoretti et al. 
2005). With BC V ~ 0.043 dLeieune k Schaerer 



20011 ). we obtain L/LQ 
2.9±0.5. 



24±8 and R/RQ 



3. Depth of the convection zone and gran- 
ulation time scales 

In Fig.[2l we compare the effective temperatures 
and radii of HD 174936 and HD 50844 with those 
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Fig. 2. — R versus T c g for metal-poor YREC mod- 
els. The radii of the bases of the surface convec- 
tion zones (grey-shaded area) are shown as dashed 
lines. 



of metal-poor evolutionary tracks, which were 
compute d with the Yale Stella r Evolution Code 
(YRE C; iGuenther et~aTl 119921 iDemarque et all 
2008) for a initial chemical composition (Y, Z) 
= (0.28, 0.01) and a mixing length parame- 
ter a = 1.8. Details ab out the model p hysics 



(2008) 



adopted can be found in iKallinger et al. 
and references therein. We estimate the masses 
of HD 174936 and HD 50844 from their positions 
in this diagram to be 1.7±0.2 and 1.75±O.25M0, 
respectively, which agr e es wi th th e findings of 
Garcia H ernandez et al. (l2009h and IPoretti et al.l 
(|2009l ) using a similar method. 

More important in this context, Fig. [2] also 
shows the radius of the base of the surface con- 
vection zone. Our intermediate-mass models have 
a very thin outer convection zone (~1 % in radius) 
during their main sequence evolution. But it starts 
to slowly expands after core hydrogen depletion in 
the H-shell burning phase before the base of the 
convection zone rapidly drops near the red bor- 
der of the instability strip and the models become 
almost fully convective on the giant branch (out- 
side the plotted range). The exact position of the 
base of the convection zone will almost certainly 
depend on the actual treatment of convection in 
the models and, in our case, on a. But the general 
picture is clear. Stars located in the 5 Scuti regime 
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of the HRD are expected to have thin, but non- 
negligible, outer convective envelopes and there- 
fore should exhibit a surface granulation signal. 
We find that the outer convection zone could be 
as deep as ~1 and 1.5% of R for HD 174936 and 
HD 50844, respectively. 

One might intuitively expect that the granu- 
lation amplitude should be very low for such thin 
convective layers. But in fact the photometric am- 
plitude of the granulation signal is not directly cor- 
related with the depth of the convection zone. It 
can be assumed that the granulation amplitude 
is inversely proportional to the total number of 
granulation cells on the stellar surface, where the 
average cell size is believed to be proportional to 
the atmospheric p ressure scale height H p (see e.g. 



Stello et alj 12007). The granulation time scales, 



on the other hand, are believed to scale with the 
ratio of the cell size and cell velocity, where the 
cells are assumed to move at speeds proportional 
to the local sound speed c s . Under the assump- 
tion of an ideal adiabatic gas, the granulation fre- 
quency (inverse time s cale), y„ rm , should scale as 



Vgran oc c s /H p (see e.g. Huber et al. 20091 and ref- 



erences therein). Unlike (5 Scuti stars, solar- type 
pulsators show very regular patterns of oscillation 
modes almost symmetrically distributed around 
the so-called "frequency of maximum oscillation 
power " , v max . According to iKjeldsen fc Bedding 
(|1995l ). v m ax is also proportional to c s /H v , and 
is believed to scale as Mj (R 2 y/T c ff). This sug- 
gests that the inverse granulation time scales de- 
termined from our global model fits scale as 



(x M ■ R- 2 ■ Tj /2 



(2) 



with the fundamental parameters in solar units. 

This relation is shown in Fig. [31 where we com- 
pare the granulation frequencies with v max of a 
sample of main-sequence and re d-giant solar-type 
pulsators ( Kallinger et al. 20091) . All parameters 
have been determined directly from CoRoT data 
using our global model approach. For the two 
5 Scuti stars we cannot of course measure v m ax 
as they do not pulsate in solar-type p-modes of 
high radial order. But we can estimate from their 
fundamental parameters what would be their v max 
if solar-type oscillations were excited in them, and 
then compare these virtual frequencies to the mea- 
sured granulation frequencies. And indeed, the 
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Fig. 3. — The characteristic v gran as a function of 
the Vmax for main-sequence and red-giant solar- 
type pulsator. Parameters on the abscissa are 
given in solar units and show how v ma x scales with 
the stellar parameters. 



two stars follow the expected linear relation. We 
therefore conclude that the granulation time scales 
are consistent with those expected from a sim- 
ple scaling relation for other stars unambiguously 
showing surface granulation, and that we have in- 
deed measured a granulation signal. 

4. Conclusions and future prospects 

We have shown that the multitude of low- 
amplitude frequencies found in the CoRoT ob- 
servations of (5 Scuti stars are consistent with 
a strongly frequency-dependent intrinsic back- 
ground signal. Such a signal is well known for 
main-sequence and red-giant stars, but has never 
been recognised in S Scuti stars. We have fit the 
power spectra of two such stars with a simple 
model, allowing us to measure the stars' granula- 
tion time scales and to disentangle the granulation 
signal from the pulsation signal. Consequently, we 
find a much smaller number of frequencies that we 
attributed to 5 Scuti-type pulsation than the hun- 
dreds to a thousand individual modes reported by 
others for the same stars. Our number of pulsa- 
tion modes is more consistent with the number of 
frequencies that can be expected from integrated 
light observations of low-degree p modes. We have 
estimated granulation time scales for the two stars 
which are consistent with those expected from a 
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simple scaling relation for stars known to exhibit 
surface granulation. 

For HD 174936, we find a low-frequency back- 
ground of ^3000 ppm 2 / /xHz, which corresponds 
to ^36 ppm in the Fourier amplitude spectrum 
of the 27-day-long CoRoT run. This should be 
roughly the same for a MOST run of comparable 
duration, and there are several S Scuti stars which 
have been monitored by MOST for this long. We 
plan to reanalyse the MOST light curves of bright 
5 Scuti stars to check for granulation signatures. 
For one year of Kepler data, we expect a signal 
of (27/365) 1 / 2 x 36 ppm ~ lOppm. For HD 50844, 
the low-frequency signal is ^125,000 ppm 2 / /iHz, 
which corresponds to ^160 ppm in the 58-day 
CoRoT run. For a 30-day MOST run, we esti- 
mate a low-frequeny background of 220 ppm, and 
for a one-year Kepler run, 60 ppm. 

Our result has shown that with the increas- 
ing quality of space-based photometry, intrinsic 
frequency- dependent noise can no longer be ne- 
glected in the frequency analyses of S Scuti stars. 
The traditional method to ascribe significances 
to candidate pulsation frequencies assuming white 
noise in a S Scuti star poses the danger of over- 
interpretating the data. This is especially impor- 
tant in the context of the upcoming Kepler ob- 
servations which will also include observations of 
8 Scuti stars, but with even higher precisions and 
longer time coverages than possible with MOST 
and CoRoT. 
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